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1. Motivation 

The CKM matrix element \Vub\ is best determined by measuring the decay rate for B — Xuiv, 
which is proportional to \Vub\^. The advantage of charmless semileptonic decays over charmless 
hadronic B decays is that the leptonic and hadronic currents of the amplitude factorize. Calculations 
of the hadronic current are difficult, since they must take into account physical mesons rather than 
free quarks. So the hadronic current is typically parameterized by form factors which can be 
calculated in the framework of QCD. 

In this analysis [jl]] \Vub\ is determined with exclusive decays which, compared with inclusive 
decays, have the advantage of reduced backgrounds but suffer from lower signal yields. We mea- 
sure A^/Aq'^, the partial branching fraction with respect to q^, the momentum-transfer squared, 
for four decay modes: B° — ;> tt^^+v, B+ — )• tt^^+v, B'^ — > p"^+v, and B+ — )- p^^+v. 

2. Data set and candidate selection 

This analysis is based on a data set of 377 million BB pairs recorded with the BaBar detector 
at the PEP-II energy-asymmetric e^e^ collider operating at the Y{4S) resonance. An additional 
sample of 35.1 fb^' of data was collected at 40 MeV below the T{4S) resonance, which is used 
to study non-BB backgrounds. Monte Carlo (MC) techniques [^] are used to simulate BB physics 
production and decay as well as detector efficiencies and resolutions. 

Signal B candidates are selected based on their three decay products: a high-momentum lepton 
(i = e^n), a hadron {n"^ ,n^\p'^ ,p^), and a neutrino. The neutrino is reconstructed from the missing 
energy and momentum in the event, and several criteria require that it be consistent with a physical 
neutrino. The second B in the event is not explicitly reconstructed; this untagged approach increases 
the signal yield but results in larger backgrounds. 

The major challenge of this analysis is to separate the relatively small signal from the much 
larger backgrounds. The largest background comes from other BB decays, in particular charmed 
semileptonic B — Xciv decays, which have a rate about 50 times that of charmless semileptonic 
B — )■ Xi,£v decays. Backgrounds also originate from 1717 {q = u,d,s,c) events, which differ from BB 
events in that they have a more jet-like topology. B — )• Xu£v decays are very similar to the signal, 
and are especially prevalent at high q^. The background composition changes as a function of q^, 
with qq contributing mostly at low and high q^, B — >• X^iv the dominant background at medium q^, 
and B — )• Xuiv the largest background at high q^. 

To reduce background contamination, neural nets are used with inputs based on event shape 
and neutrino kinematics. Nets are trained separately in different ranges of q^ against each of the 
three dominant backgrounds: qq, B — )• Xc£v, and B — )• X,/v. 

Detailed control sample studies are used to check the agreement between simulation and data. 
The data and MC agree well in S — )• Xc£v-enhanced samples and exclusively reconstructed B^ — )• 
[)*-£+ V decays, suggesting that the neutrino reconstruction is well-modeled in the simulation. 

3. Signal yield fit and A^/Aq^ 

The signal yield is extracted with an extended binned maximum likelihood fit [Q] that takes into 
account the statistical uncertainties of both the data and MC samples. The fit is performed to the 
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three-dimensional distribution of AE, mes, and q^. The variables AE and mns test the consistency 
of the candidate with a B decay, and are defined as: AE = Eg — \/s/2 and wies = \J{y/^jTp^^-p^ , 
where ^/s is the center-of-mass energy of the colliding beams, and Eg and p*^ are the center-of-mass 
energy and momentum of the recontructed B. The momentum-transfer squared, q^, is calculated 
as the square of the sum of the lepton and neutrino 4-momenta. To improve the resolution, the 
neutrino momentum is scaled to set AE to zero in the calculation of q^. 

The four modes are fit simultaneously with isospin constraints for the charged and neutral 
modes. The fit yields 10604 ± 376 signal B n£v decays and 3332 ± 286 signal B plv decays, 
from which the branching fractions are calculated: 

=^(5° ^ TT^^+v) = (1.41 ± 0.05 ± 0.07) X 10""^ , 
^{B^^p-tv) = (1.75 ±0.15 ±0.27) x 10"^ , 

where the first error is statistical and the second error is systematic. Separate measurements of the 
branching fractions from single-mode fits to charged or neutral B samples are found to be consistent 
within statistical uncertainties with the result from the combined four-mode fit. 

The error on each of these branching fractions is dominated by systematic uncertainties. For 
B — )■ Tiiv, the largest systematic uncertainties arise from the spectrum of Kl, which carry away 
missing energy and momentum and thus impact the neutrino resolution; and from the shape of the 
fit distributions for the qq background. For B — )• p£v, the largest uncertainties arise from the shape 
function parameters and branching ratio of non-resonant B X^iv decays, which cannot be easily 
distinguished from the signal. 




Figure 1: Partial branching fraction of B*^ — > tt £^v (left), ^ p t^v (right), in bins of Values 
measured from data (points with errors) are compared with predictions from light cone sum rules (LCSRl 
[Bp and LCSR2 [§ for B^ tt^^+V, and LCSR [0] for ^ p"^+v), unquenched lattice QCD (HPQCD 
[ ^ I), and quark-model (ISGW2 [^) theory calculations, normalized to equal area. The BGL parameterization 
[ ^ is shown with the shaded region indicating the uncertainty on the fitted parameters. For B^ — > n^t^v, 
the dashed lines indicate an extrapolation of the prediction to the full <7^ range. 

The B — )• %lv signal yield is fit in 6 bins of c^, and the B ^ plv yield is measured in 3 bins 
of as shown in Figure after correction for efficiency, detector resolution, bremsstrahlung. 
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and radiative effects. The measured A^/Aq^ distribution is compared with several theoretical 
form-factor calculations. For — )• n^i^v the agreement with the data is best for the HPQCD 
lattice calculation [||], but the lattice QCD predictions are only valid for > 16GeV^, the earlier 
light cone sum rules (LCSR) calculation (LCSR 1 [|]) for q'^ < 16GeV^, and the more recent LCSR 
calculation (LCSR 2 [^]) for q^ < 12GeV^. For — )• p^£^v, the branching fraction measurements 
are not precise enough to discriminate between the predictions from LCSR and the ISGW2 
quark-model calculation [^. 

4. Determination of I I 

The value of \Vub\ can be determined from the measured partial branching fraction over a 
limited range of q^ and an integral of the form factor over the same range of q^: 

Vuh\ - \\ ^ 2 ' AL,[q^i^,q^^) - — j / Pj,\J+[q j dq , 

where Tq = (1.530 ±0.009) ps is the B° lifetime. 

\Vub\ is determined from form-factor theoretical predictions from light cone sum rules at low 
q^ and lattice QCD at high q^, and in each case the theoretical uncertainty of about 15% dominates. 
These measurements of \Vub\ are shown in the first three lines of Table |T| 





q^ Range 


A.^ 


A? 


\Vuh\ 




(GeV^) 


(10-4) 


(ps-i) 


(10-3) 


LCSR 1 [|] 


0-16 


1.10 ±0.07 


5.44±1.43 


3.63 ±0.121°:^^ 


LCSR 2 [|] 


0-12 


0.88 ±0.06 


4 00+' °i 


3.78±0.13l°-^5 


HPQCD [U 


16-26.4 


0.32 ±0.03 


2.02±0.55 


3.21±0.17+°:^5 


FNAL/MILC [|ll|] 


0-26.4 


1.41 ±0.09 




2.95 ±0.31 


HPQCD [§ 


- 26.4 


1.41 ±0.09 




2.99 ±0.35 



Table 1: |V„/,| derived from B %lv decays for various regions and form-factor calculations. The first 
three values use data points from only a limited range, while the final two values use data from the full 
range. Experimental and theoretical errors on | V„f,| are listed separately for partial-g'^ determinations and 
combined into a single error for the fuU-g'^ determinations. 

\yuh\ can also be determined using a combined fit to data and theoretical predictions over the 
full range of . Three fit parameters describe a quadratic polynomial in the BGL form-factor 
ansatz [|l0|], and a fourth parameter defines the relative normalization of theory and data, which is 
proportional to \yuh\- The advantage of this method is that it reduces the theoretical uncertainty. The 
relative error contributions are 3% from the branching-fraction measurement, 5% from the shape of 
the spectrum determined from data, and 8.5% from the form-factor normalization obtained from 
theory. The result of this fit is shown in Figure ^ and the values of | V„6 1 for two lattice form-factor 
calculations are shown in the last two lines of Table Unfortunately, the lattice calculations are 
only valid at high (f, where the decay rate is lowest and the experimental uncertainties are largest. 
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In summary, \Vub \ has been determined with two methods. Using data from a limited range, 
larger values of \Vub\ are found with theory predictions for low than with predictions for high 
q^. A combined fit to theory and data from the full q^ range reduces the overall error by a factor of 
two. The value of \Vub\ from both methods is smaller than most determinations of \Viih\ based on 
inclusive B — )- Xu£v decays, which are in the range (4.0 — 4.5) x 10^^. 




q2 (GeV^) 



Figure 2: Partial branching fraction of Z?" n^i'^v in bins of q^, with the combined fit to data (black 
points) and lattice theory (squares or triangles) represented by the solid red line. Not all theory points are 
used because of large correlations. The shaded region indicates the uncertainty on the fitted parameters. 
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